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Abstract 

We examine electrically and color neutral quark matter in /^-equilibrium focusing 
on the possibility of mixed phases between different color superconducting phases. 
To that end we apply the Gibbs criterion to ensure phase equilibrium and discuss the 
external conditions under which these mixed phases can occur. Neglecting surface 
and Coulomb effects we find a rich structure of different mixed phases with up to four 
components, including 2SC and CFL matter as well as more "exotic" components, 
like a phase with us- and <is-pairing but without uti-pairing. Preliminary estimates 
indicate, however, that the mixed phases become unstable if surface and Coulomb 
effects are included. 



1 Introduction 

The possible existence of color superconducting matter in the core of a neutron star is 
subject to intensive discussions since it has been discovered § that at intermediate 
densities the energy gaps can be of the order of 100 MeV. These large values have been 
obtained in the scalar color antitriplet channel where instanton interactions || as well 
as interactions derived from single-gluon exchange predict the largest attraction. There 
are two different condensation patterns in this channel, depending on whether or not the 
strange quarks, which are more massive than the light up and down quarks, participate 
in forming a condensate, 

SAA> = C75 T A X A , 1p) . (1) 

Here both, ta and \a> are the antisymmetric generators of SU(3), i.e., the antisymmetric 
Gell-Mann matrices (A, A' e {2,5,7}), acting in flavor and color space, respectively. In 
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the two-flavor color superconducting phase (2SC) where only the light quarks are involved 
in the condensation, the flavor index in Eq. (|l|) is restricted to A = 2. In this case it is 
always possible, without loss of generality, to perform a color rotation such that the 2SC 
phase is described by s 2 2 ^ and saa> = if (A, A') ^ (2,2). The color-flavor-locked 
(CFL) phase [|J is characterized by nonvanishing condensates S22, S55, and S77, i.e., it 
contains ud as well as us and ds pairs. In the limiting case of equal masses and chemical 
potentials for all quarks these three condensates are equal, whereas in general they can 
be different from each other. 

Most of the calculations so far have been performed using a common chemical potential 
for all quarks (see, e.g., |5], U). It has been found that in that case a large region of the 
phase diagram is occupied by a 2SC phase |7|, §]. For the interior of neutron stars, it is 
however important, to consider electrically and color neutralQ matter in /3-equilibrium. 
As pointed out by Alford and Rajagopal this strongly disfavors or even rules out the 2SC 
phase in compact stars [H|. For example consider a system of massless up and down 



quarks together with electrons, but without strange quarks. Since the density of electrons 
is small (see, e.g., [[H]]), to achieve electric neutrality the density of (i-quarks must be 
almost twice as large as the density of u-quarks, and hence /i^ ~ 2 1 / 3 /i u . This means that, 
e.g., for fi u = 400 MeV, the Fermi momenta of u and d differ by about 100 MeV, making 
ud BCS-pairing very difficult. The presence of negatively charged strange quarks therefore 
is likely to favor pairing even if the strange quark mass is considerably larger. Alford and 
Rajagopal tackled that problem, performing an expansion in terms of the strange quark 
mass. They found that, whenever the 2SC phase is more favored than no pairing at all, 
the CFL phase is even more favored. 



Recently this has been reinvestigated by Steiner, Reddy and Prakash [12], who ana- 



lyzed electrically and color neutral quark matter within an NJL-type model, taking into 
account density dependent quark masses self-consistently. For those conditions which are 
relevant for neutron stars older than a few minutes, i.e., cold matter without trapped 
neutrinos, the authors found the 2SC phase to be most favored within a small region 
of quark number chemical potentials. They claim, however, that this region is likely to 
disappear, once the hadronic phase is properly taken into account [jEJ. 

The analyses of Refs. |L0], [12| are based on a comparison of homogeneous phases: 
For a given quark number chemical potential /i the authors search for the phase which 
maximizes the pressure after the constraints of electrical and color neutrality have been 
imposed. The latter is achieved by introducing an additional chemical potential /iq which 
is related to electric charge and two chemical potentials /i 3 and /i 8 which are related to 
color. The values of these new chemical potentials needed to neutralize matter for a fixed 
value of /i depend on the phase. For instance the normal phase, where color invariance is 
not spontaneously broken, is color neutral for ^3 = fig = 0, whereas a nonvanishing value 
of \iq is needed for electric neutrality. The situation is quite opposite to the CFL phase 



where /z 8 7^ and fj,Q — pi 3 — are needed fl2| . Finally, in the 2SC phase both, /i 8 and 



/iq, must be nonzero. Hence neutral matter in one phase is never in chemical equilibrium 

*Strictly speaking, color neutrality is not sufficient, but color singletness has to be imposed. This 
does, however, not induce a large cost in energy such that we can consider matter which is only color 
neutral. 
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with neutral matter in a different phase, even if their quark number chemical potentials 
fx are the same. In particular the points of equal pressure do not fulfill the Gibbs criteria 
for a phase transition stating that the pressure and all chemical potentials must be the 
same in coexisting phases. 

One possibility to resolve this problem is to give up the requirement of separately 
neutral phases and to consider mixed phases of two components in chemical equilibrium 
which are only neutral in total. This procedure has been pushed forward by Glendenning 
in the context of the quark-hadron phase transition in neutron stars where a similar 



problem occurs ||13|| . Since Glendenning did not consider color-superconducting phases, he 
only had to care about electric neutrality. In this case a neutral mixed phase can obviously 
be constructed in those regions of the phase boundary where the charge densities of the 
two components have opposite signs. 

In the present paper we generalize this procedure to construct electrically and color 
neutral mixed phases. This is more difficult because for two charged and colored com- 
ponents a mixture which is electrically neutral is in general not color neutral and vice 
versa. In the four- dimensional space spanned by lm, /Iq, /z 3 and /i 8 the phase bound- 
aries are three-dimensional hypersurfaces. As we will show the regions where electrically 
and color neutral mixed phases are possible correspond to one- dimensional lines on these 
hypersurfaces. 

Our calculations are performed within the same NJL type model which was the basis of 
the analysis of Ref. [EI. It includes interactions which allow for the condensation patterns 
described above as well as for quark-antiquark condensates (ipip). The latter, describing 
dynamical chiral symmetry breaking in the vacuum, in particular generate dynamically 
a constituent strange quark mass which can be of the order of the quark chemical po- 
tential We restrict ourselves to bulk quark matter in mean-field approximation. Of 
course if mixed phases are formed the domains belonging to a single component must have 
finite sizes determined by the electric and color-electric fields and the surface tension. In 
fact, if the latter is too large the mixed phase will become unstable and a single charged 
interface between the two phases will be more favorite. This scenario has been suggested 
for the transition from the hadronic phase to the CFL phase WM. A crystalline color 



superconducting phase |15| in the transition region would be another possibility. 

The paper is organized as follows: In Sec. ^| we introduce the conserved charges in 
our system and the related chemical potentials. We discuss the Gibbs condition for phase 
boundaries. Thereafter, in Sec. |||, we present the model which is used to describe quark 
matter. We continue in Sec. ^ with an examination of the resulting phase diagram in 
the space of the various chemical potentials. In Sec. |5| we construct electrically and 
color neutral mixed phases and discuss the resulting phase structure of neutral quark 
matter. Finally, in Sec. |6| we summarize and discuss some open questions related to our 
investigations. 
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2 Conserved charges, chemical potentials and phase 
transitions 



We consider a system of quarks, electrons and muons with the number densities n/ )C , n e , 
and n M , respectively. Here / = u,d,s ("up", "down", "strange") refers to the flavor and 
c = r,g,b ("red", "green", "blue") refers to the color. The total flavor and color densities 
are then given by 

/ 

As mentioned in the Introduction we are mainly interested in describing the conditions 
present in compact stars older than a few minutes, when neutrinos can freely leave the sys- 
tem. In this case lepton number is not conserved and we have four independent conserved 
charges, namely the total electric charge 

2 11 

nQ = -n u - -n d - -n s ~n e -n^ (3) 
and the three color charges n r , n„, and rib or, equivalently, their linear combinations 



n = n r + rig + n b , n 3 = n r - n g , n 8 = -j={n r + n g — 2n b ) . (4) 

V3 



Here n corresponds to the total quark number density, while % and n s describe color 
asymmetries. Note that n/3 also describes the conserved baryon number. The four 
conserved charges are related to four independent chemical potentials, /i, /x 3 , /x 8 , and \xq. 
The individual chemical potentials /i/ jC "felt" by the quarks of flavor / and color c are 
then given by corresponding components of the diagonal matrix 

1 1 

\L = ii + n Q (-r 3 + ^=ts) + /i 3 A 3 + ^ 8 A 8 . (5) 

Here, as before, and Xj are operators in flavor space and color space, respectively. The 
electron and muon chemical potentials are simply fi e = /i M = — fiQ- This implies 

Hd,c = Hs,c = Hu,c + He for all c , (6) 

which is usually referred to as /^-equilibrium. 

The total thermodynamic potential of our model is given by the sum of a quark part 
and a leptonic part, 

fi(T,{/i,}; X ) = n q {T,{Hi}; X ) + ^i{T,Hq) , (7) 

where {/ij} = {/i, yU 3 , /x 8 , hq}- For the leptonic part we simply take a gas of noninteracting 
massless electrons and massive muons. The quark part will be further specified in Sec. |3|. 
Here we only note, that it depends on a set of condensates which we denote by x- m 
general there are different solutions for x a t given temperature and chemical potentials. 
Then the stable one is the solution which corresponds to the lowest value of Q. At 



4 



the boundaries which separate two different phases the corresponding thermodynamic 
potentials coincide, 

n(T,{fi t }; X i) = fi(T,{/i l }; X 2) , (8) 

where Xi 7^ X2 if the phase transition is first order. 

In a given phase the various densities defined above can be obtained from the thermo- 
dynamic potential as 

on 

We are mostly interested in electrically and color neutral matter, which is characterized 

by 

n Q =n 3 = n 8 = . (10) 

For homogeneous phases the consequences of these constraints have been discussed in 
Refs. [jTUl O. The goal of the present article is to analyze possible mixed phases which 



could exist along first-order phase boundaries, Eq. (||). 

Suppose we have a mixed phase consisting of two components, 1 and 2. In general 
the densities n\ and resulting from Eq. (0) will not be identical in both phases. In 



particular, in general Eq. fllO| ) will not be fulfilled simultaneously for both components. 
However, as indicated in the Introduction it is sufficient to demand that the average 
charge and color densities of the mixed phase vanish Jl3| . If the two components occupy 



the volume fractions x\ and x% — 1 — xi, respectively, the average densities are given by 

rii = xinf' + (1 — X\)nf' . (11) 



This is zero for 



(2) 

Xl = "72) (IT ' ( 12 ) 
n) — n) 



To be meaningful the solution must be in the interval < x± < 1. This is fulfilled when 
the charge densities n[ and nf* have opposite signs, which is an obvious prerequisite for 
a charge neutral mixture. For a single charge, e.g., uq, it is the only one. However, in 
order to get simultaneous neutrality for three charges, Eq. (|T0|) , we have to require that 
the result of Eq. ( |I2"D is the same for % — Q, 3, and 8. This is the case when 

(1) (1) (1) (2) (2) (2) / 1Q x 

In our numerical calculations we will restrict ourselves to T = 0, relevant for compact 
stars older than a few minutes. Then the phase boundaries, Eq. (§), are three-dimensional 
surfaces in the four-dimensional space spanned by fi, /ig, /i3, and [1$. Since Eq. ([131 ) 
imposes two additional constraints, electrically and color neutral mixed phases can be 
constructed along a one- dimensional line. In the simplest case this line starts at a point 
where the neutrality line of phase 1 (n^ = = = 0), meets the phase boundary 
and it ends where the neutrality line of phase 2 meets the phase boundary. Between these 
two points Xi changes continuously from 1 to 0. However, if the neutrality line meets 
another phase boundary before, one has three coexisting phases, 

Q(T,{fH};xi) = tt(T,{in};X2) = V{T, {^}; * 3 ) . (14) 
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In this case the neutrality condition reads 



/ (1) (2) ( 

/ nW nW fir 



Nx = 



n 



\ 



n 



n 



n 



(1) „(2) 
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(!) „(2) „(3) 
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f *1 \ 











(15) 



In order to find a non-trivial solution for x, we must have det N = 0. Together with 
Eq. ([14]), this again restricts the possible solutions to a one- dimensional subspace. More- 
over, since the fractions X{ should not be negative, for each m = Q, 3, 8 the densities nf 
must not have the same sign for all % = 1, 2, 3. 

Finally, there could even be a mixed phase, consisting of four components. The cor- 
responding phase boundary is one- dimensional and again the region of possible neutral 
mixed phases is further restricted by the requirement that the various fractions x-i should 
not be negative. 



3 Quark model 

The quark part of our model is defined by the Lagrangian 

^e// = 4>(i$-m)?P + Cgg + Cgg, (16) 

where ip denotes a quark field with three flavors and three colors. The mass matrix 
rh has the form rh = diag(m u ,md,Tn B ) in flavor space. Throughout this paper we will 
assume m u = rrid, whereas the strange quark mass m s will be different. To study the 
interplay between the color-superconducting diquark condensates and the quark- ant iquark 
condensates we consider an NJL-type interaction with a quark-quark part 

= H £ ^t^T A \ A ,C^ T )^ T Cz l5 T A \ A ,^). (17) 

A=2,5,7 A' =2,5,7 

and a quark-antiquark part 

8 

As = G^2[(M) 2 + (^75^) 2 ] - tffdet/^l+Tsty) + det/ (^(1-78)^ 

a=0 



Here r = y 1 1/ is proportional to the unit matrix in flavor space. Eq. fll8|) corresponds 
to a typical 3-flavor NJL-model Lagrangian which has often been used to study meson 
spectra [|l^, [HJ or properties of quark matter at finite densities or temperatures |18|, [19], [11 



It consists of a ?7(3)x, x ?7(3)^-symmetric 4-point interaction and a 't Hooft-type 6-point 
interaction which breaks the the U A (1) symmetry. The above interaction terms might arise 
from some underlying more microscopic theory and are understood to be used at mean- 
field level in Hartree approximation. In particular we do not consider any contribution 
from the 6-point interaction to the diquark condensate. The same model has also been 
employed in Ref. |12[ and for the case of one common chemical potential in Ref. ||. 



6 



In order to calculate the mean-field thermodynamic potential fi g at temperature T 
and the various chemical potentials we linearize £ e ff m the presence of the three quark- 
antiquark condensates 0/ = (//), / = u,d,s, and the three diquark condensates saa, 
A = 2,5, 7. Introducing the constituent quark masses 



Mi = mi - AG(j>i + 2K(fi j( fi k 
and the diquark gaps 



k) = any permutation of (u, d, s) 



(19) 



A A = -2Hs AA , (20) 
and employing Nambu-Gorkov-formalism the result can be written in the following way: 

d 3 p 



+ 2G(£ 



(2tt) 3 

+ cfi 2 d + 4> 2 S ) - 4K<fi u , 



+ H(\s 22 \ 2 + \s 55 \ 2 + \ S77 \ 



(21) 



Here x — {4>u, (fid, 0s, S22, S55, s 77 } denotes the set of condensates. S 1 is the inverse 
fermion propagator 



s-\p) 



j> - M + £7° A 2 7 5 r 2 A2 + A 5 75r 5 A5 + A 7 j 5 t 7 \ 7 

-A* 75 r 2 A 2 - A;7 5 r 5 A 5 - A; 75 r 7 A 7 j> - M - /x 7 ° 

(22) 

where M = diag(M u , Md, M s ). In Eq. (|2~ID , it has to be evaluated at p — (iu n ,p), where 
u> n = (2n — 1)tvT are fermionic Matsubara frequencies. 

For the case of equal chemical potentials an explicit expression for Q q has been dis- 
cussed in detail in Ref. 0. For unequal chemical potentials this is more involved and we 
only state the general structure as far as necessary for our later discussion. As shown in 
Refs. [fTO , 12] , pOl S^ 1 can be decomposed into several blocks. Six of them have a 2 x 2 
structure in the combined color-flavor space, 



M h 



C1J2C2 



M 



h 



A75 



-A* 7 5 f f 



M f 



(23) 



J hfi2 - lvl h 

where p^ c = ft ± /i/, c 7o- These blocks describe pairing of two species of quarks with 
flavors fi and f 2 and colors C\ and c 2 , respectively. More precisely, one gets the following 
three combinations [F2|: First, green up quarks are paired with red down quarks, where 
A = A 2 , second, blue down quarks are paired with green strange quarks (A = A 5 ), and 
third, red strange quarks are paired with blue up quarks (A = A7). The other three 2x2 
blocks can be generated simply by interchanging the role of particles and antiparticles 
and lead to the same combinations. 

The remaining contribution to Q has the following 6 x 6-structure: 



M 



uds 



( ~ M u 











-A275 


-A575 \ 





i>X 9 - M d 





-A275 





-A 7 7 5 








ft - M s 


-A 5 7 5 


-A 7 7 5 








A^7s 
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Kr ~ Mu 
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A 7 7s 





K 9 ~ Md 





V A* 75 
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K b -m s ) 



(24) 
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involving red w-quarks, green <i-quarks and blue s-quarks. This represents the particle 



and antiparticle contributions to the 3 x 3-block mentioned in Ref. |20 . 

As demonstrated by Rajagopal and Wilczek for a simplified two-quark model, the 
pairing of two quark species forces their densities to be equal |24]]. Originally, it was 



expected that in the CFL phase, where all quarks are paired, all densities n^ c are equal. 
This would mean that the CFL phase is always electrically and color neutral. However, 
as pointed out by Steiner, Reddy, and Prakash , only those quarks which are paired in 
the same 2x2 block have the same density, whereas the densities could differ for different 
blocks. Furthermore the argument does not apply to the 6x6 block. For the CFL phase 
this means n U)9 = rid t r, n Ut b = n Str , and n^b = n s ^ g , and thus 

n u = n r , n d = n g , n s = n b (CFL) . (25) 

This relation guarantees neutrality of CFL matter under the rotated electromagnetism 



Q |]T0], 11211 ) but in general it does not preclude the presence of ordinary electric or color 
charges. Note, however, that color neutral CFL matter is automatically electrically neutral 
as long as no leptons are present. 

In the 2SC phase, where A 5 = A 7 = 0, Eq. (|24D can be further decomposed, and we 
obtain a new 2x2 block involving red w-quarks and green <i-quarks. Together with the 
other 2x2 block which contains A 2 this leads to the relations 

n u ,r = n d>g , n U:9 = n djr (2SC) . (26) 

The corresponding relations for other phases, e.g., with two nonvanishing diquark con- 
densates, can be obtained analogously. 

The further elaboration of the thermodynamic potential contains only straight forward 
manipulations, which will not be presented here. In Eq. ( ^Tl) Q q depends on x, i.e., on 
the choice of the condensates 0/ and saa- In order to find the self-consistent solutions we 
have to look for the stationary points of the potential with respect to these condensates. 
As mentioned before, there is often more than one selfconsistent solution. In this case the 
stable solution is the one which minimizes the thermodynamic potential. 

Our model is incomplete without specifying the parameters to be used in the numerical 
calculations. To this end we adopt the parameters of Ref. |19| , which have been determined 



by fitting the pseudoscalar meson masses and decay constants in vacuum: a sharp 3- 
momentum cutoff A = 602.3 MeV to regularize the integrals, the bare quark masses 
m u = rrid = 5.5 MeV and m s = 140.7 MeV, and the coupling constants G = 1.835A -2 
and K = 12.36A~ 5 . Similarly, one could fix the quark-quark coupling constant H by 
fitting baryon masses within a Fadeev approach |2l], For simplicity, however, we 
follow Ref. |T2[ (preprint versions vl and v2) and take H = G.Q This is similar in 
magnitude to the value we have chosen in Ref. @ where we took H = 3/4 G, but with a 
larger Gf\ 



'la the latest version of Ref. |12| (v3), the authors take H = 3/4 G. Qualitatively, this does not change 
their results. Therefore, in the following we will not distinguish between the different versions. However, 
for any quantitative comparison with our results, the reader should refer to version v2. 

Hn Ref. M we used a four-point interaction with the quantum numbers of a single-gluon exchange, 



s 




300 400 500 300 400 500 

H [MeV] |i [MeV] 

Figure 1: Gap parameters and densities at T = ^ = //§ = fj,Q = as functions of \i. Left: 
M u = M d (solid), M s (dashed), A 2 (dotted), and A 5 = A 7 (dashed-dotted). Right: n/10 (solid), 
n§ (dashed), and tiq (dotted). 

4 Phase structure 

In order to construct the equation of state for electrically and color neutral quark matter 
(including possible mixed phases) we first have to explore the phase structure of the model 
and the charge densities as functions of the chemical potentials /ij. 

We begin with the standard case of a single, color and flavor independent, chemical 
potential // for all quarks, i.e., /i 3 = /j 8 = fiq = 0. For the parameters specified above 
we obtain the results which are displayed in Fig. |IJ. In the left panel the constituent 
masses and diquark gaps are shown. One can clearly distinguish three phases separated 
by first-order phase boundaries: At low chemical potentials we find a normal phase with 
vanishing diquark condensates, followed by the 2SC phase with A 2 ^ and finally the 
CFL phase where also A5 and A7 are nonzero. At the phase boundaries we observe strong 
discontinuities in the quark masses ]7[ 

In the right panel we show the corresponding densities. Note that the quark number 
density n (solid line) has been divided by 10 to fit to the scale. The dotted line corresponds 
to the electric charge density Uq, the dashed line to the color density n s . The color density 
n 3 is identically zero. As one can see, all densities vanish in the "normal phase", i.e., this 
phase corresponds to the vacuum. In fact, it has to be like this, because as soon as up and 
down quarks are present, their Fermi surfaces are subject to a Cooper instability leading 
to the formation of the diquark condensate s 22 . (This argument will no longer go through, 
once we have switched on one of the other chemical potentials which lift the degeneracy 

without six-point interaction. In this case a simple Fierz rearrangement yields H : G — 3 : 4. This 
relation also holds for a two-flavor instanton interaction. Hence, in our present model a simple estimate 
for the value of the diquark interaction, e.g., in the wd-channel, is given by H — 3/4 G e //, where G e ff = 
G — K(ss). In vacuum this is slightly larger than G, but it is a density dependent quantity and differs 
for different phases. Such effects have not been taken into account in our present work, but they are 
certainly worthwhile to be studied in future. 
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of the Fermi surfaces of all up and down quarks.) 

The two other phases carry both, electric and color charges. The electric charge of the 
2SC phase is easily understood. Since /iq = 0, there are no leptons and the densities of up 
and down quarks are equal. Moreover, in our example there are no strange quarks, which 
are too heavy to be populated in this regime. Hence the total electric charge density is 
given by nq = n/6. The nonvanishing color density n 8 reflects the fact that for equal 
chemical potentials the densities of the paired (red and green) quarks are larger than the 
density of the unpaired (blue) quarks [23]. Numerically we find (n r — n&)/n = 10% at the 
lower boundary and (n r — = 3% at the upper boundary of the 2SC phase. 

Just above the transition to the CFL phase this ratio does not change very much, 
whereas the electric charge density drops significantly due to a strong increase of the 
density of strange quarks. To a large extent, this is caused by a sudden drop of the 
strange quark mass, but this is only part of the story. For instance, at /i = 500 MeV 
we have M u = = 17.2 MeV and M s = 216.5 MeV. Using these numbers in a free 
gas approximation we would expect nq = 0.049 fm -3 , whereas numerically we find nq = 
0.030 fm -3 . This difference is caused by the diquark pairing, which links the flavor 
densities in the CFL phase directly to the color densities, as discussed in Eq. (|25|). For 
n 3 = one finds nq = l/(2y / 3)^8; i n agreement with our numerical results. 

Electrically and color neutral CFL and 2SC matter has been constructed in Ref. ||12|| . 
According to the charge densities n 8 and nq which have to be removed, this in general 
requires the introduction of nonvanishing chemical potentials /ig and fiq. In order to see 
how these additional chemical potentials influence the phase boundaries, we first study 
them separately. 

In Fig. ^ we show the phase diagram in the /i — /ig-plane for fig = /i 3 = 0. The 
(first-order) phase boundaries are indicated by solid lines. We find again the three phases 
discussed before, i.e., the normal phase, the 2SC phase, and the CFL phase. For /x 8 = 
we have seen that the "normal phase" actually corresponds to the vacuum. However, 
when fib = fj, — 2/v / 3/is becomes larger than the vacuum masses of the light quarks (the 
region above the dashed line), blue up and down quark states can be populated forming a 
gas of unpaired blue quarks. Here we have neglected that in principle these quarks could 
pair in a different channel 0, |26|. In any case we should keep in mind that our model 
is not suited for a realistic description of the low-density regime, where confinement and 
hadronic degrees of freedom have to be taken into account. 

For the color superconducting phases we have also indicated the lines of color neutral 
matter (dotted). As discussed in the previous section, in the CFL phase color neutral 
quark matter is automatically electrically neutral as well, i.e., in the CFL phase the dotted 



line corresponds to the neutral matter solution, which has been determined in Ref. fl2 |. 
It meets the phase boundary to the 2SC phase at \i = 465.7 MeV and /i 8 = —32.5 MeV. 
The 2SC matter which is in chemical and mechanical equilibrium with the neutral CFL 
matter at this point carries both, electric and color charge, nq = 0.464 fm -3 and n s = 
—0.329 fm -3 . In the next section, this point will be the end point of the 2SC-CFL mixed 
phase. Unlike color neutral CFL matter, color neutral 2SC matter is not electrically 
neutral but positively charged. In fact, a nonvanishing /ig does not change the ratio of 
up and down quarks and hence, as long as no strange quarks are present, nq/n = 1/6 as 
before. 
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Figure 2: Phase diagram in the u — /ig plane for T = f/,3 = fJ,Q = 0. The various phases separated 
by the solid lines are characterized by different non-vanishing diquark gaps A, as indicated in 
the figure. In the non-superconducting phase quarks are present only above the dashed line. 
The dotted lines are the lines of color neutral matter. In the CFL phase this also corresponds 
to electrically neutral matter. 



In Fig. |3] we show the phase diagram in the \i — yUg-plane for /x 8 = /i 3 = 0. Since 
we are interested in neutralizing the electrically positive 2SC phase, we choose \iq to be 
negative. As long as this is not too large, we find again the normal phase at lower values 
of /1, the 2SC phase in the intermediate region and the CFL phase for large fi. This 
changes dramatically around \xq ~ 180 MeV where both, the 2SC phase and the CFL 
phase disappear and a new phase emerges. This phase is analogous to the 2SC phase but 
with ds pairing, instead of ud pairing ("2SCd s "). In a small intermediate regime there is 
even another phase which contains us and ds but no ud pairs ("SC us+ d s "). 

Qualitatively, the existence of these phases is quite plausible: At low values of \/iq\ the 
Fermi momenta of the up and down quarks are relatively similar to each other, whereas 
the strange quarks are suppressed because of their larger mass. With increasing negative 
/iq, however, the up quarks become more and more disfavored and eventually the Fermi 
momenta are ordered as p F < p s F < p F . It is then easy to imagine that only ds pairing or 
- in some intermediate regime - us and ds pairing is possible. 

Following this argument one might expect that there is always a value of /j,q, where 
the Fermi momenta of up and strange quarks are equal and hence the 2SC phase should 
be either followed by the CFL phase or by a phase with its-pairing only. However, this is 
not the case because of the discontinuous behavior of the quark masses. This is illustrated 
in Fig. |] where the diquark gaps and constituent quark masses are shown as functions of 
\iq for fixed \i = 390 MeV and /13 = fi$ = 0. The 2SC-SC ns+ d s phase transition takes 
place at hq = —178.6 MeV, corresponding to /x u = /x + 2/3/ig ~ 270 MeV and = 
fi s = fi — 1/3/iQ ~ 450 MeV. Below the transition point the strange quark mass is even 
larger than 460 MeV and, consequently, no strange quarks are present. At the transition 
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Figure 3: Phase diagram in the fJL—fiQ plane for T = = u$ = 0. The various phases separated 
by the solid lines are characterized by different non-vanishing diquark gaps A, as indicated in 
the figure. In the non-superconducting phase quarks are present only above the dashed line. The 
"+" and "-" signs indicate the sign of the electric charge density in the corresponding region. 
The dotted line corresponds to electrically (but not color) neutral matter in the CFL phase. 



point the strange quark mass drops to 310 MeV and the nominal Fermi momentum 
p s F = — Ml is immediately larger than p F . 



The stability of the various condensates is rather well described by the criterion |24] 



Pf ~ Pf 



V2 



A*- , (27) 



where A^- is the gap related to the pairing of the quark species i and j. In the 2SC 
phase just below the phase boundary we have A2 = 132.8 MeV, slightly larger than A 2 = 
127.4 MeV. At the phase boundary A 2 rises to 133.6 MeV due to a sudden increase of the 
up quark mass by more than 40 MeV. Taking the earlier value of A 2 , the above criterion 
is no longer fulfilled, which is consistent with our numerical result that the w<i-pairs break 
up. This level of agreement is certainly better than one should expect. In fact, in the 
SC us+ ds phase we find A 5 continuously decreasing from 50.8 MeV to 49.1 MeV whereas 
A5 increases from 48.2 MeV to 52.6 MeV, slightly violating Eq. (|27|) . It is nevertheless 
possible to understand the break-up of the us pairs, which occurs at \iq = —183.0 MeV, 
from the fact that at this point A5 jumps to 62.6 MeV due to a further increase of M u 
and a further decrease of M s . Moreover, the fact that we always find A5 « A5, at least 
in this example, indicates that the SC us+ d s phase is rather fragile and might disappear 
upon small variations of the model parameters. 

In the phase diagram, Fig. [| we also indicated the sign of the electric charge density 
for the various regions, and the line of electrically neutral matter in the CFL phase (dotted 
line). Note that there is no other electrically neutral regime in this diagram (apart from 
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Figure 4: Diquark gaps and quark masses for T = ^3 = fJLs = 0, and [i = 390 MeV as functions 
of fiQ. Left: A2 (solid), A5 (dashed), and A7 (dotted). Right: M u (solid), (dashed), and 
M s (dotted). 



the vacuum at small /1 and [iq = 0). In the normal phase, there are again no quarks below 
the dashed line, corresponding to the line ji — 1/3/xq = This region is nevertheless 
negatively charged due to the leptons which are present for any /ig < 0. Above the dashed 
line there are also d quarks rendering the matter even more negative. (In the right corner 
of this phase there is also a very small fraction of u quarks.) 

The "new" phases, 2SCd s and SC us+ rf s , are also negatively charged. On the contrary, 
the entire 2SC phase is positively charged, even at the largest values of \(jlq\. This is 
illustrated in Fig. where the various charge densities n.j divided by the total quark 
number density n are plotted as functions of fiQ, again for fixed /i = 390 MeV and 
A*3 — A*8 = 0. As expected, uq/u (solid line) decreases with increasing negative Hq. 
However, in the 2SC phase (0 > \lq > —178.6 MeV) it stays positive and before the point 
of neutrality is reached the phase transition to the SC us+( fe phase takes place. 

The difficulty to obtain electrically neutral 2SC matter can be traced back to the fact 
that according to Eq. (|26D , the sum of red and green u quarks is equal to the sum of red 
and green d quarks. As long as no strange quarks are present, the related positive net 
charge can only be compensated by the blue quarks and the leptons, which would require 
a very large negative [iq. However, before this point is reached it becomes more favorite 
to form a different phase with a relatively large fraction of strange quarks which then also 



phase 


N 


2SC 


2SC US 


2SCd s 








CFL 


diquark gaps 




A 2 


A 5 


A 7 


A 2 , A 5 


A 2 , A 7 


A 5 , A 7 


A 2 , A 5 , A 7 



Table 1: Phases and corresponding non- vanishing diquark gaps. 
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0, and \x = 390 MeV as functions of [j,q: hq/h 



participate in a diquark condensate.0 As discussed above, the selfconsistent treatment, 
which leads to a sudden drop of the strange quark mass and hence to a sudden increase 
of the strange Fermi momentum, is crucial in this context. 

In the above discussion we have not considered the effect of a non-vanishing chemical 
potential //§ on top of a non- vanishing fig. Since the blue quarks are the main carriers 
of negative electric charge in the 2SC phase, one could hope that increasing the number 
of blue quarks, as necessary for color neutrality, could also help to electrically neutralize 
2SC matter. It turns out, however, that the rather small values of fis which are needed 
for color neutrality (see Ref. (T2| and Fig. [2J) do not change the above results qualitatively. 

So far we always kept yU 3 = 0. In fact, there is no need to vary /i 3 , as long as we 
are mainly interested in finding electrically and color neutral solutions of homogeneous 
normal, 2SC, or CFL matter. However, as we will see in the next section, the construction 
of neutral mixed phases requires also non- vanishing values of /i 3 . In this context we will 
encounter another phase, which is not present in Figs. and [| For illustration we consider 
a plane in the four-dimensional {/Zj}-space where \x and fiQ are taken as independent 
variables and /i 3 and /i 8 are given by /i 3 = —Hq/2 and /i 8 = — /xq/7 — 30 MeV. The 
relevance of this particular choice will become more clear in the next section. Here we 
just note that /z 3 = — /xq/2 means that fi u>r = fid, g • Also the sum ^ s ,r+^u,b, corresponding 
to the chemical potential related to a pair of a red strange quark and a blue up quark, 
equals the sum // S)S + fM^b, corresponding to the chemical potential related to a pair of 
green strange quarks and blue down quarks. Together with the relations given above 
Eq. (p5|) and the isospin symmetry of the original Lagrangian this implies for the CFL 
phase that n u = n d or, equivalently, n r = n g and thus n 3 = 0. 

In Fig. |6| we show a small part of the resulting phase diagram. Here, in addition to 
the standard 2SC and CFL phases, we find a phase where only u and s quarks are paired 

§This is very similar to the arguments of Alford and Rajagopal [To| . The main difference is that we 
do not compare different neutral phases with each other, but phases in chemical equilibrium. 
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Figure 6: Phase diagram in a plane denned by /i and fiQ as independent variables and 
^3 = — hq/2 and ^8 = — /Uq/7 — 30 MeV. The various phases separated by the solid lines 
are characterized by different non-vanishing diquark gaps Aj as indicated in the figure. 



("2SC US "). This means, we have found examples for all three A^, A = 2,5,7, being the 
only non-vanishing scalar diquark gaps in some regime. Taking all possible combinations 
of no, one, two, or three of these condensates (see Table HD, the phases SC u d +us and 
SC u d+ds, be., the combinations A 2 + A 5 and A 2 + A 7 are the only ones we have not 
encountered. These phases might exist as well, but we did not perform a systematic 
investigation of this question, which was not the goal of this work. 



5 Mixed phases 

Having gained some overview about the general phase structure of our system we are 
now in the position to construct electrically and color neutral mixed phases. Our starting 
point is ii — 465.7 MeV, /is = —32.5 MeV, and ^ — fj,g — where the line of neutral 
CFL matter meets the boundary to the 2SC phase. At lower values of /i, mixed phases 
become possible and are energetically favored as long as Coulomb and surface effects are 
neglected. Following the guidelines worked out in Sec. || we find nine regimes characterized 
by different compositions of coexisting phases (see Table [2J). The corresponding chemical 
potentials /ij as functions of /i are displayed in Fig. [7|. In Table || we also list the corre- 
sponding minimal and maximal quark number densities, averaged over the components 
of the respective mixed phase. 

In the regime closest to the region of homogeneous neutral CFL matter (430.6 MeV < 
fi < 465.7 MeV), we find a mixed phase consisting of a CFL component and a 2SC 
component. The volume fraction X2sc of the 2SC component is displayed in the left panel 
of Fig. |8|. In the higher-/^ part of this region it is completely negligible, but even at the 
lower end it remains below 2%. Consequently, the CFL component must stay almost 



15 



neutral by itself. Indeed, the relative charge densities rii/n, i = 3, 8, Q, (right panel of 
Fig. H) are very small. As we have discussed in the previous section, n 3 -neutrality of the 
CFL phase is maintained by the relation ^3 = — /aq/2. For the actual values of ^3 and 
/iq/2 in the 2SC-CFL mixed phase we find a deviation of less than 1% from this relation, 
while /is can approximately be fitted by /is = — A*q/7 — 30 MeV. This is the reason why 
we have calculated the phase diagram shown in Fig. ^ with these constraints. 

As can be seen there, the 2SC-CFL phase boundary meets the boundary to the 2SC US 
phase at \x = 430.6 MeV and we get a three-component neutral mixed phase, consisting of 
2SC, CFL and 2SC MS . Below that, on a short interval in /1, we even find a four-component 
neutral mixed phase (2SC, CFL, 2SC US , and SC us+ d s ) before upon further decreasing /t 
the system goes over into a neutral 2SC-CFL-SC us+ d s mixed phase. 

In Fig. |9] the volume fractions of the various components of the mixed phases are 
plotted as functions of /i. Whereas the 2SC-CFL-mixed phase (Fig. ||) is completely dom- 
inated by the CFL component, thereafter the CFL fraction becomes quickly smaller with 
decreasing /1, while in particular the 2SC component, becomes more and more important. 
At jj, = 407.7 MeV the CFL component disappears completely. 

An admixture of normal quark matter is found below /i = 388.8 MeV. The fractions 
of the superconducting phases other than the 2SC phase then rapidly become smaller 
and vanish at fi = 388.6 MeV, while the fraction of normal matter strongly increases. [] 
Nevertheless the 2SC phase stays the dominant component for /a > 360 MeV. 

As discussed in the previous section, apart from the vacuum there is no region of 
stable neutral non-superconducting quark matter. Therefore the normal-2SC mixed phase 
cannot end in normal homogeneous quark matter but only in the vacuum. To that end, the 
chemical potential \xq must finally go to zero. Eventually, at /i = 348.6 MeV, |/iq| drops 
below 60 MeV and we enter the regime where the normal phase only consists of electrons 
without quarks (cf. Fig. |3|). This means the corresponding mixed phase consists of 
(electrically positive) 2SC-droplets surrounded by regions without quarks and neutralized 

^Note, however, that all volume fractions Xi are continous functions in the entire mixed phase region. 



components 


fi [MeV] 


n [fm 3 ] 


N, 2SC 


340.9 - 388.6 


0.00 - 1.56 


N, 2SC, SC MS _|_<2 S 


388.6 - 388.7 


1.56 - 1.56 


N, 2SC, SC us +ds> 2SC US 


388.7 - 388.8 


1.56 - 1.56 


2SC, SC us +ds, 2SC ns 


388.8 - 395.4 


1.56 - 1.74 


2SC, SC us +ds 


395.4 - 407.7 


1.74 - 2.59 


2SC, SC us +ds, CFL 


407.7 - 426.5 


2.59 - 2.90 


2SC, SC us +dsi CFL, 2SC MS 


426.5 - 427.1 


2.90 - 2.93 


2SC, CFL, 2SC US 


427.1 - 430.6 


2.93 - 3.11 


2SC, CFL 


430.6 - 465.7 


3.11 - 3.92 



Table 2: Composition of electrically and color neutral mixed phases, corresponding quark 
number chemical potentials and average quark number densities. The various components are 
defined in Table |l|. 



16 



200 



150 

> 

S 100 

50 





350 



375 



400 

|^ [MeV] 



425 



450 



Figure 7: Chemical potentials ^ corresponding to the electrically and color neutral mixed 
phases, listed in Table ||: —[Iq (solid), ^3 (dashed), and — (dotted). 
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Figure 8: Quantities related to the neutral 2SC-CFL mixed phase as functions of the quark 
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Figure 9: Volume fractions of the various components in the mixed phase region as functions 
of fi: normal (thin solid), 2SC (thick solid), CFL (dashed-dotted), SC us+ d s (dashed), 2SC US 
(dotted). 



by a homogeneous background of electrons. Since the electrons are color neutral, the 2SC 
component must be color neutral by itself. This is maintained by an increase of |/xg| in 
this regime. At fi — 340.9 MeV we finally reach the vacuum. 



Qualitatively our results are in agreement with those obtained in Ref. [12|, where no 



mixed phase was admitted. The regions which are dominated by the CFL phase and 
the 2SC phase, respectively, correspond roughly to those values of the chemical potential 
where the authors of Ref. jT2|| find that the CFL or the 2SC phase, respectively, is fa- 
vored. Of course, since we started from the same Lagrangian, these homogeneous neutral 
solutions also exist in our model. We find, for instance, below \i = 465.7 MeV a color 
and electrically neutral CFL phase which is, however, less favored than the 2SC-CFL 
mixed phase. At this point we should recall that in our calculations we have neglected 
the surface energy and the energies of the electric and color-electric fields. Clearly the 
mixed phases are only stable if the gain in bulk free energy is larger than these so-far 



neglected contributions. If this is not the case we recover the results of Ref. [|12]. In a 
neutron star we would then find homogeneous phases separated by sharp interfaces. 

For the two-component mixed phases the surface and Coulomb contributions can be 
estimated along the lines of the analysis performed in Ref. flM for the interface between 
nuclear and CFL matter. Because of the smallness of the 2SC fraction in the 2SC-CFL 
mixed phase the gain in bulk free energy compared with a pure CFL phase is extremely 
small in this case (less than 0.1 MeV/fm 3 ). It is therefore more or less obvious that 
this mixed phase will not be stable. The situation looks a little more promising for the 
2SC-SC us+( 2 S -phase at fi « 400 MeV where the rarer phase, i.e., the SC us+( fe-phase, is 
more populated. Here we find a gain in bulk energy of about 4 MeV/fm 3 . However, even 
in this case a calculation along the lines of Ref. shows that this gain is weight out 
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by Coulomb and surface energy already for relatively small values of the surface tension 
er pa 10 MeV/fm 2 . Note that the true value of a is an unknown quantity which does not 
follow from our model. However, as argued in Ref. [14] it is likely to be much larger than 
the above value. 

For mixed phases with three or four components the situation is obviously more com- 
plicated. It seems, however, rather unlikely that an admixture of more than two compo- 
nents can be favored if the surface and Coloumb energy already inhibits the existence of 
a phase with two components. The more components we have the more contributions to 
surface, Coulomb and color charge energy we have. On the other hand, even if, e.g., a 
3-component mixed phase turns out to be unstable, we cannot exclude to find a stable 
2-component mixed phase instead. A definite answer to these questions lies beyond the 
scope of this paper and we postpone it to future work. 



6 Summary and Discussion 

We performed an investigation of quark matter under conditions which possibly occur 
in the interior of neutron stars older than a few minutes, i.e., dense, electrically and 
color neutral matter in /3-equilibrium. Recently it was shown by Alford and Rajagopal 
that these constraints have important effects, in particular that they might disfavor the 
standard two-flavor color superconducting phase (2SC) in neutron stars [jnj. In the 
present article we focussed on the formation of mixed phases, i.e., different coexisting 
phases in chemical equilibrium. In this case the neutrality conditions have only to be 
fulfilled globally @. 

Our model calculations are performed within an NJL-type model, which allows to 
treat diquark condensates on an equal footing with quark- ant iquark condensates, leading 
to density (and phase-) dependent effective quark masses f7j. In accordance with our ex- 
pectations the regimes of lower or higher quark number chemical potentials are dominated 
by normal (non-superconducting) matter and by matter in the color-flavor locked (CFL) 
phase, respectively. In addition, we find a certain regime of quark number chemical po- 
tentials where the 2SC phase is the dominant component. Naively one would expect that 
the above three phases represent the only constituents of the mixed phases we encounter. 
However, in addition we find admixtures of more exotic phases, like, e.g., a SC us+ d s -phase 
or a 2SC us -phase. These phases occur because the chemical potentials, in particular fiq, 
necessary to enforce electric neutrality, are often large enough such that the ordering of 
the Fermi momenta is reversed: p F <p s F < p F . 

Simple estimates show that the gain in bulk energy due to the formation of these mixed 
phases is probably not large enough to compensate for the energy arising from surface, 
and electric and color charge effects. This seems to indicate that homogeneous phases of 
neutral quark matter separated by sharp interfaces are favored. This could include the 
possibility of a stable neutral 2SC phase. However, a more careful analysis is certainly 
needed before a definite answer to this question can be given. There are other aspects 
which deserve to be studied more closely in future: 

At low densities we certainly need a more realistic description of the hadronic phase, 
which in the present model at best is represented by a vacuum or quark phase with 
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spontaneously broken chiral symmetry. A more complete calculation, including a realistic 
equation of state of nuclear matter, has to show whether in particular the N-2SC-phase 
which we find at rather low densities is robust. 

At higher densities we should also take into account the possibility of a CFL-phase 
with kaon condensation p9| , |30| . Recently the authors of Ref. Jl2| have estimated the 



effect of this possibility on the equation of state of homogeneous neutral quark matter. 
They found that the kaon condensate lowers the critical quark chemical potential for the 
2SC-CFL phase transition by about 16 MeV. This is not dramatic but could make a 
difference if for instance the transition point to the hadronic phase turns out to be in the 
same region. 

We also neglected the possibility of the formation of crystalline (LOFF) phases [|15|, [27 



These become eventually favored over phases with BCS pairing if the Fermi surfaces 
are pulled apart. As we are mostly dealing with pairing of quarks with unequal Fermi 
momenta, we suppose that especially in the vicinity of phase boundaries, the possible 
formation of a crystalline phase should be taken into account. In fact, the occurrence 
of mixed phases which are necessarily related to non-uniform structures in space already 
points into this direction. The neglect of gradient terms in our mean-field approach is only 
justified if the structures are much larger than the average distance between the quarks. 
Unfortunately, we cannot say much about the sizes in the mixed phases as long as we 
do not know the surface tension. However, taking relatively small surface tensions as 
required to keep the mixed phases energetically stable, the sizes of the rarer components 
are typically only a few times larger than the inter-quark distances. 

Recently another interesting pairing configuration has been discussed by Liu and 
Wilczek ||28|| . A so-called "interior gap" can be formed in systems with unequal Fermi 



momenta where the Fermi momentum of the heavier fermion species is larger than that 



of the lighter species. The authors of Ref. claim that the "interior gap" configuration 
is favored if 

p a F + p b F M a 
1 > 2p F M a + M b ' (28) 

where M a < M b and p F < p b F . Usually this possibility can be excluded in studies of 
superconducting quark matter, since the Fermi momentum of the heavier strange quarks 
is smaller than the Fermi momenta of up and down quarks. But, as we have discussed 
above, the existence of "exotic" phases, like, e.g., the SC us+( fe-phase, can be traced back 
to the fact that the Fermi momenta are ordered as p F <p s F < p d F . Therefore, us- 
pairing in the "interior gap" configuration seems possible. Let us look at some numerical 
example: At (J, « 400 MeV we find SC us+ds matter with M u » 60 MeV, M s « 300 MeV, 
p F ~ 285 MeV, and p s F ~ 350 MeV. This leads to values of about 0.2 for the right hand 
side of Eq. Thus it seems probable that such a phase occurs. 

It has also to be mentioned that we have restricted ourselves to diquark condensates 
saa' with A = A'. Originally this was motivated by the fact that, e.g., in the 2SC phase 
any linear combination of S22, S25, and S27 can be rotated into the 22-direction by a global 
color transformation. However, in principle we can always construct a color neutral mixed 
phase without applying color chemical potentials by combining several components of the 
same condensate but rotated into different color directions. Within our present formalism 
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this would further lower the free energy, but it is rather unlikely, whether this is still 
the case when the color electric energy of the various domains are taken into account. 
Anyway, problems, like the difficulty to find stable 2SC matter, mostly arise from the 
constraint of electric neutrality, whereas imposing color neutrality is only a minor effect. 

Note that we only considered the situation relevant for neutron stars older than a 
few minutes, i.e., when the neutrinos can freely leave the star. The authors of Ref. [jP2]] 
performed an extended analysis looking also at a system at finite temperature and lepton 
content. These complications are likely to be relevant during the evolution of the star from 
a proto-neutron star at high temperature where neutrinos are trapped to a cold compact 
star where neutrinos can freely leave the star. Steiner, Reddy and Prakash flllj showed 



that in this case the 2SC phase is favored because the CFL phase excludes electrons and 
can therefore not easily accommodate a finite lepton number chemical potential. A study 
of mixed phases in this context is left to future work. 
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